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1 Introduction 

o 



HIPPOPO (Hadron Induced Particle Production Off Pomeron Outpour) is a program 
to simulate the production of hadrons via the multipomeron exchange mechanism in a 
high energy hadron collision. The underlying algorithm is based upon the Quark Gluon 
O ' String Model developed over the years by many authors 0-[H. 



Or 



o 



c3 



For the initial colliding beams the User can choose among the following species: nu- 
cleons, atomic nuclei, tt, K, cf) and J/ip mesons, A, £ and H hyperons (of any explicitly 
stated charge), and their antiparticles. Also, the program allows photon beams. In the 
latter case, the photon is treated as a superposition of p, and J / ip mesons according to 
' the Vector Dominance Model. 

<3\ 



The produced hadron may be a tt, r\, p, uj, K, <p, D, D s or J/ifj meson, nucleon N, A, 
S, H or Q hyperon, A c or S c baryon, an antiparticle or an excited state of any mentioned 
hadron. 

The program yields the fully differential inclusive cross section. However, it is not a 
full event generator, because it only focuses on the inclusive production properties of the 



user-defined hadron, while the accompanying particles are not taken into account. 

2 Physics input 



X 

An interaction between the colliding hadrons is assumed to proceed by means of the 
multipomeron exchange mechanism. The particular contribution to the cross section due 
to the exchange of n Pomerons reads [Qj: 

'.«>= £(l-«-g£). W 

with 

£ = In s, a P = 8vt7p e ?A , z = 2C 7p [R 2 + a^]" 1 s A , 

and C, 7p, R 2 , a' P and A being the model parameters (see section 3.3 for their numer- 
ical values). Also, there exist contributions from the elastic and diffraction dissociation 
processes, a e i and u DD , which correspond to the absense of Pomerons (n — 0): 

*o(0 = *p[f(f )-/(*)] with f(z) = T,—^^ ^ 
a el = {I/O) *„(£), ° DD = (1 " ^o(0- (3) 
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Summing up, the total inelastic cross section is: 

oo 

Vtot = V DD + J2 (T r l - ( 4 ) 



n=l 



According to the Quark Gluon String Model, each Pomeron is treated as a pair of 
colour strings attached to the partons in the colliding hadrons. The fragmentation of 
colour strings results in the production of particles. The inclusive differential cross section 
for a hadron of type h reads: 

where the particle production from Pomeron strings is described by the string fragmen- 
tation functions <p n : 

ipt = a h [FUx + ,n)F^(x.,n) + F^(x + ,n)F v h al (x-,n) 

+(n - 1) [Fg ea (x+, n)F^x.,n) + F^(x + , n)F* ea {x„,n)\], (6) 

with x± = \[{x\ + x 2 F ) l l 2 ±xf], x± = 2m\l 'y/s, so that x + X- = (m^) 2 /s and 

X-\- — X— — Xp. 

Each term in (|6|) corresponds to an independent string stretched between oppositely 
coloured partons from different beams. The contribution from each string is given by 
a product of two independent factors, Fi(x + ) and Fj(x-), related to the endpoint partons 
(moving in the positive and negative directions, respectively). In turn, the factors F(x±) 
are given by the convolution of universal functions: 

F h (x ± , n) = J2 t f*( x '> n ) G ?( x ±/ X '> Pt) T(x f /x', Pt , n) dx>, (7) 
i Jx ± 

where fi(x, n) stands for the probability to find a parton of type i carrying the momentum 
fraction x in the beam hadron, and Gf(x,p±) stands for the probability that the fragmen- 
tation of a type i parton yields a hadron h with the longitudinal momentum fraction x 
and transverse momentum p±. The parton distribution functions fi(x,n) are specific for 
a given beam particle, and the fragmentation functions G'l(x,p±) are specific for a given 
outstate hadron. 

The diffraction dissociation term was taken in the form 



<p* D = ^F\x + , 1) + y/*IF h (x-, 1) (8) 
if the fragmenting and the produced particles are both baryons or both antibaryons, and 

= (3/2)[ v ^T + y/xj F h (x + , 1) F h (x^ 1) (9) 



DD 

in all other cases. 

The parton distribution functions have the (Reggeon theory inspired) generic form: 

f(x,n)=x a t(l-x) fi t+ n - 1 (10) 

with the exponents «/ and f3f determined by the flavour content of the beam particle. 
Similarly, the generic form of the parton fragmentation functions is: 

G(x,p T ) = x a °{l - x f G+2a '^T (11) 
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with the exponents «g and (3q determined by the flavour content of the produced hadron. 
The parameters a/, /3f, o>g and fie are expressed in terms of the relevant Regge intercepts 
and represent the basic model assumptions. At the same time, to provide an interpolation 
between the limiting points x — > and x — > 1, some extra factors may be introduced in the 
functions G^(x,pt)- The latter ones are only of phenomenological meaning and cannot 
be derived from 'first principles'. The full list of parton distribution and fragmentation 
functions is too long to be presented in this paper, but these functions are easily readable 
from the FORTRAN code. 

Finally, the weight function T(z,pj<,n) in (0) is present for the hadron transverse 
momentum distribution. According to |13| , it is parametrized in the form: 

T(z,p T ,n) = [b 2 z /2n(l + b z m h )] exp{-b z (m h ± - m h )] (12) 

with b z = 7V [1 + (2 — ^/ n )pz 2 ], and 7^ and p being still new phenomenological parame- 
ters. 

3 Program components 

3.1 General structure and external references 

The program structure is divided into several pieces. These are the user job cards 
and the initiating routines (collected in the file run.f); the table of parton distribution 
functions (file parton . f ); the table of fragmentation functions (files f ragml . f , f ragm2 . f , 
fragm3.f); the particle production algorithm as described in Sect. 2 (file model. f); the 
Monte Carlo integration routine VEGAS [^TJ with its intrinsic random number generator 
(file vegas . f ). 

When evaluating the convolution in ([7]) the program refers to the standard function 
DGAUSS from CERN library. The normalizing factors for parton distribution functions are 
calculated with the use of the DGAMMA function, also from CERN library. 

For histogramming purposes, the program uses the standard PAW package. 

3.2 Subroutines and Functions 

SUBROUTINE CKNAME ( i 1 , i2 , in) 

Checks the consistency of the user given names for the beam and the produced particles. 
SUBROUTINE PARAM 

Introduces the model parameters. 
SUBROUTINE HYBRID 

Assignes model dependent intercepts to hybrid and excited states. 
SUBROUTINE MIX 

An algorithm to include composite beams (Atomic nuclei, for example) and hadrons 
of indefinite quark content (short and long living Kaons, 7/ mesons, etc.). 

SUBROUTINE CXNPOM(sqs ,nmax) 

Calculates the multipomeron cross sections, equ. ([I]). 
FUNCTION F0(z) 

Returns the f(z) value, equ. (|2|). 



3 



FUNCTION FXN(xveg,wgt) 

The integrand expression for VEGAS []2l|, with xveg being the array of phase space 
variables and wgt the weight factor due to the particular phase space binning (supplied 
by VEGAS). 

FUNCTION Fldif(x), F2dif(x) 

1st and 2nd beam diffractive dissociation functions. 
FUNCTION Flval(x), F2val(x) 

1st and 2nd beam valent colour endpoint functions. 
FUNCTION Flbar(x), F2bar(x) 

1st and 2nd beam valent anticolour endpoint functions. 
FUNCTION Flv(x) , F2v(x) 

1st and 2nd beam sea colour endpoint functions. 
FUNCTION Flb(x), F2b(x) 

1st and 2nd beam sea anticolour endpoint functions. 
FUNCTION T(zF,n,hmass,pt2) 

Transverse momentum distribution function. 



SUBROUTINE NORM(nmax) 

Normalizes all parton distribution functions. 
FUNCTION Bfun(ex,el) 

Euler's Beta function. 



FUNCTION 

FvalU (x , n , beam) , FvalD (x , n , beam) , FvalS (x , n , beam) , FvalC (x , n , beam) 

Valent u, d, s, c quark distribution functions. 
FUNCTION 

FbarU (x , n , beam) , FbarD (x,n, beam) , FbarS (x , n , beam) , FbarC (x , n , beam) 

Valent u, d, s, c antiquark distribution functions. 
FUNCTION 

FseaU(x,n,beam) , FseaD(x,n,beam) , FseaS(x,n,beam) , FseaC(x,n,beam) 

Sea u, d, s, c quark distribution functions, apply to antiquarks also. 
FUNCTION 

FvalUU(x,n,beam) , FvalUD(x,n,beam) , FvalDD(x,n,beam) 

Valent uu, ud, dd diquark distribution functions. 
FUNCTION 

FvalUS(x,n,beam) , FvalDS(x,n,beam) , FvalSS(x,n,beam) 

Valent us, ds, ss diquark distribution functions. 
FUNCTION 

FbarUU(x,n,beam) , FbarUD(x,n,beam) , FbarDD(x,n,beam) 

Valent uu, ud, dd antidiquark distribution functions. 
FUNCTION 

FbarUS (x , n , beam) , FbarDS (x , n , beam) , FbarSS (x , n , beam) 

Valent us, ds, ss antidiquark distribution functions. In the above functions: 

x is the momentum fraction, 

n is the number of exchanged Pomerons, 

beam is the beam hadron the distribution functions refer to. 
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FUNCTION 

DvalU(x,hadr,pt2) , DvalD(x,hadr ,pt2) , DvalS (x , hadr ,pt2) , DvalC(x,hadr ,pt2) 

The u, d, s, c quark diffraction dissociation functions. 
FUNCTION 

DbarU(x,hadr,pt2) , DbarD(x,hadr ,pt2) , DbarS(x,hadr ,pt2) , DbarC(x,hadr ,pt2) 

The u, d, s, c antiquark diffraction dissociation functions. 
FUNCTION 

DvalUU(x,hadr,pt2) , DvalUD(x,hadr,pt2) , DvalDD (x , hadr ,pt2) 

The mm, ud, dd diquark diffraction dissociation functions. 
FUNCTION 

DvalUS(x, hadr.pt 2) , DvalDS(x,hadr,pt2) , DvalSS(x,hadr,pt2) 

The us, ds, ss diquark diffraction dissociation functions. 
FUNCTION 

DbarUU (x , hadr , pt2) , DbarUD(x,hadr,pt2) , DbarDD(x,hadr,pt2) 

The mm, ud, dd antidiquark diffraction dissociation functions. 
FUNCTION 

DbarUS^,hadr_ L pt2) , DbarDS(x,hadr,pt2) , DbarSS(x,hadr,pt2) 

The us, ds, ss antidiquark diffraction dissociation functions. 

FUNCTION 

GvalU(x,hadr,pt2) , GvalD(x,hadr ,pt2) , GvalS(x,hadr ,pt2) , GvalC(x,hadr ,pt2) 

The m, d, s, c quark fragmentation functions. 
FUNCTION 

GbarU(x,hadr,pt2) , GbarD(x,hadr ,pt2) , GbarS(x,hadr ,pt2) , GbarC(x,hadr,pt2) 

The m, d, s, c antiquark fragmentation functions. 
FUNCTION 

GvalUU(x,hadr,pt2) , GvalUD(x,hadr,pt2) , GvalDD (x , hadr ,pt2) 

The mm, ud, dd diquark fragmentation functions. 
FUNCTION 

GvalUS(x,hadr,pt2) , GvalDS(x,hadr ,pt2) , GvalSS(x,hadr,pt2) 

The ms, ds, ss diquark fragmentation functions. 
FUNCTION 

GbarUU(x_,hadr L pt2) , GbarUD(x,hadr,pt2) , GbarDD (x , hadr ,pt2) 

The mm, ud, dd antidiquark fragmentation functions. 
FUNCTION 

GbarUS^,hadr_ L pt2) , GbarDS(x,hadr,pt2) , GbarSS(x,hadr,pt2) 

The ms, ds, ss antidiquark fragmentation functions. In the above functions: 

x is the momentum fraction, 

hadr is the hadron to appear from fragmentation, 

pt2 is its transverse momentum squared. 

SUBROUTINE DDMODE 

Establishes the incoming and outcoming particle types as 'meson', 'baryon' or 'an- 
tibaryon', that is important for diffraction dissociation contributions. 

SUBROUTINE ABSRAT 

Establishes the absolute production rate. 
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SUBROUTINE VEGAS (FXN , AVGI , SD , CHI2A) 

Perform multidimensional Monte Carlo integration 
DOUBLE PRECISION FUNCTION RANDOM (SEED) 

Random number generator W2 . 



SUBROUTINE WRIOUT(hwgt) 

A user supplied routine to manage output information, hwgt is the correctly normalized 
weight that reproduces the physical cross section (|5]) when is summed over all the 
generated events. 



3.3 Common blocks 



C0MM0N/INPUT/beaml,beam2, hadron (All are CHARACTER* 6 names) 

beaml 1st beam particle type (defined by the User). 

beam2 2nd beam particle type (defined by the User). 

hadron The produced hadron type (defined by the User). 
C0MM0N/BEAMS/beaml,beam2, hadron (All are CHARACTER* 6 names) 

beaml 1st beam particle type (for internall use). 

beam2 2nd beam particle type (for internall use). 

hadron The produced hadron type (for internall use). 
C0MM0N/AT0MN/Atoml(2) ,Atom2(2) 

Atoml(l) Atomic number of the 1st beam (if nucleus). 

Atoml(2) Atomic weight of the 1st beam (if nucleus). 

Atom2(l) Atomic number of the 2nd beam (if nucleus). 

Atom2(2) Atomic weight of the 2nd beam (if nucleus). 
COMMON/EXCIT/Level 

The outgoing hadron excitation level, Level=0 is the ground state. 

C0MM0N/VDMES/M1 (3) ,M2(3) 

The keys that switch on and off the vector meson component in a photon. Ml(i) and 

M2(i) refer to the 1st and the 2nd beams, respectively, and i=l,2,3 stand for p, <fi 

and J/ip meson contributions. 

COMMON/KINEM/sqs , xl , x2 , xF , y , hmass , pt2 



SC L S Total c.m.s. energy, ^/s [GeV]. 

xl Parton momentum fraction in the 1st beam, x + 

x2 Parton momentum fraction in the 2nd beam, X- 

xF Feynman variable of the produced hadron, xf- 

y Rapidity of the produced hadron, y. 

hmass Mass of the produced hadron, m# [GeV]. 

pt2 Transverse momentum squared, p\ [GeV 2 ]. 



COMMON/STRNG/n , nmax 

The number of Pomerons exchanged and the maximal number of Pomerons 
COMMON/POMER/Coef f , Gammp , Rpom2 , AlPom , Delta 



Coeff 


= 1 


5 


C = l + (r DD / (J eh equ.© 


Gammp 


= 3 


64 


Pomeron residue parameter, 7p [GeV -2 


Rpom2 


= 3 


56 


Pomeron size parameter, R 2 [GeV -2 ] 


AlPom 


= 


25 


Pomeron trajectory slope, a' P [GeV -2 ] 


Delta 


= 


07 


Pomeron intercept criticality, A 
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The above parameters are fixed in the model and are not recommended to be altered. 



COMMON/CXPOM/CXN (0 : NP) 

The multipomeron exchange cross sections <J DD , u n . 

COMMON/REGGE/alp , aR , aP , a J , aN , aD , aV , aX , aO , aE 

The intercept parameters for Reggeon, p, ip, J /if), N, A, A, E, Q trajectories. Please 
note that the parameters aN , aD , aV , aX , aO must not be literally copied from Par- 
ticle Data. They do not have actual meaning of baryon intercepts, but only enter the 
model in artificial combinations that correspond to fictituous trajectories of multiquark 
hybrid states. aE is the effective intercept shift for excited states (a model). 

COMMON/TRNSV/gamh , gamr 

Transverse momentum distribution parameters 7 h and p, see equ. (12). 

COMMON/STRAN/Snucl , Cnucl 

Strange and Charm sea suppression parameters. 

COMMON/PIPHI/Cpi , Cph , Cps 

Fragmentation parameters for ir and <fi mesons. Note that only the associated pro- 
duction of the type <pKK is considered for <p mesons. The OZI-violating nonplanar 

diagrams are not included at present. 
COMMON/PNLAM/Cnl , Cn2 , Cll , C12 

Fragmentation parameters for p, n and A. 
COMMON/KAONS/Ck , Cks , Ckl , Ck2 

Fragmentation parameters for Kaons. 
COMMON/HYPER/Csl , Cs2 , Cxi , Cx2 , Col , Co2 

Fragmentation parameters for E, S and Q hyperons. 
COMMON/CHARM/Cd , Cdl , Cd2 , Ccl , Cc2 , Cf , Cf s , Cf c , C j , C j c 

Fragmentation parameters for D, D s , J/if> mesons and A c baryon. Note that only 

the associated production of the type J/ipDD is considered for J/if> mesons. The 

OZI-violating nonplanar diagrams are not included at present. 
COMMON/SUPER/Cxc , Cxc 1 , Cxc2 

Fragmentation parameters for S c baryons. 
COMMON/VDPAR/VDM (3) 

Photon to vector meson coupling constannts, Aira/fy. 
COMMON/RATES/rrho , reta , retap , romm , rexcit 
COMMON/YIELD/ Ctotal 

Normalization parameters for hadron production rates. 

C0MM0N/INDEX/ex,el 

The exponents that parameterize parton distributions. 
C0MM0N/CN0RM/CvalU(2,NP) , CvalD(2,NP), CvalS(2,NP), CvalC(2,NP) 

,CbarU(2,NP) , CbarD(2,NP), CbarS(2,NP), CbarC(2,NP) 

,CseaU(2,NP) , CseaD(2,NP), CseaS(2,NP), CseaC(2,NP) 

,CvalUU(2,NP) ,CvalUD(2,NP) ,CvalDD(2,NP) 

, CvalUS (2 , NP) , CvalDS (2 , NP) , CvalSS (2 , NP) 

,CbarUU(2,NP) ,CbarUD(2,NP) ,CbarDD(2,NP) 

, CbarUS (2 , NP) , CbarDS (2 , NP) , CbarSS (2 , NP) 
The normalization factors for parton distributions in the 1st and the 2nd beams, for 

an n-Pomeron exchange. 



7 



C0MM0N/BVEG1/XL(10) ,XU(10) ,ACC 

XL (10) Array of lower limits for phase space variables 

XU(10) Array of upper limits for phase space variables 

ACC Accuracy parameter 

COMMON/BVEGG/NDIM , NCALL , ITMX , NPRN 

NDIM Number of phase space dimensions 

NCALL Number of random points per iteration 

ITMX Maximal number of iterations 

NPRN Print /noprint parameter 

C0MM0N/SEED/NUM 

NUM Random number 

3.4 Job cards 

To run the program, the User has first to initiate the random number generator by 
the card NUM=1, to set the total c.m.s. energy and to specify the colliding beams and the 
hadron to be produced. 

Also, the User has to define the lower XL(i) and the upper XU(i) limits for the 
independent variables that parametrize the phase space, i.e. lnpf, and y. (Note the use 
of logarithm of the transverse momentum.) The lower limit of the transverse momentum 
is an arbitrary small number (it should be only nonzero, to avoid formal arithmetical 
conflicts). The upper limit should be chosen in a reasonable agreement with the total 
c.m.s. energy sqs. 

The wanted histograms have also to be booked in the MAIN program. 
It is recommended to perform calculations in two steps. A short preliminary run 
optimizes the VEGAS grid to the integrand function shape: 

NCALL = 1000 ! number of points per iteration 

ITMX = 5 ! number of iterations 

NPRN =0 ! do not fill histograms 

CALL VEGAS (FXN , AVGI , SD , CHI2A) 
After that one can start a long run to accumulate large statistics: 

NCALL = 200000 ! number of points per iteration 

ITMX = 1 ! number of iterations 

NPRN =1 ! do fill histograms 

CALL VEGAS1 (FXN, AVGI, SD,CHI2A) 
The quantity that is to be plotted in histograms for the physical cross section is given by 
hwgt (see SUBROUTINE WRI0UT). 
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